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ABSTRACT 

We present the first analysis of an optical to mid-infrared spectrum of the T7.5 dwarf Gliese 570D 
with model atmospheres, synthetic spectra, and brown dwarf evolution sequences. We obtain precise 
values for the basic parameters of Gl 570D: T e e = 800 - 820K, logg(cm/s 2 ) = 5.09 - 5.23, and 
logL/ 'Lq = —5.525 to —5.551. The Spitzer IRS spectrum shows prominent features of ammonia (NH3) 
that can only be fitted by reducing the abundance of NH3 by about one order of magnitude from the 
value obtained with chemical equilibrium models. We model departures from chemical equilibrium 
in the atmosphere of Gl 570D by considering the kinetics of nitrogen and carbon chemistry in the 
presence of vertical mixing. The resulting model spectrum reproduces the data very well. 
Subject headings: stars: abundances — stars: atmospheres — stars: individual (Gliese 570D) — stars: 
low-mass, brown dwarfs 
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1. INTRODUCTION 

Ammonia is the only significant nitrogen-bearing com- 
pound that can be readily observed in brown dwarfs. Its 
presence in very cool brown dwarf atmos p heres has been 
expec ted for some time |Tsuii| 119641 iM^riej^e^aLl 
19961 IFeglev fc Lodders \ 199 6t ISaumon et all 
2000l lAhard et all 120011 lLodders k, Feglev I 120021 
Burrows Sudarskv. fc Limine! l200 ^ ISaumon et al. I 
2003). ISaumon et all l)2000ft reported a detection of 
weak NH3 features in the near infrared spectrum of the 
T7p dwarf Gl 229B but a clear detection had to await 
the Spitzer InfraRed Spectrograph (IRS) observatio n 
of the T1/T6 binary e Indi Bab ijRoellig et al. 1 120041. 
IRS spectra also reveal the presence of NH3 in all 
observed dwarfs with spectral types of T2 and later 
l|Cushing et al. ll2006Th 

In addition to its importance as a strong absorber 
in very cool brown dwarfs, NH3 provides an important 
window into the phenomenon of non-equilibrium chem- 
istry. Because some chemical reactions proceed very 
slowly at the relatively low temperatures encountered 
in the atmospheres of T dwarfs, vertical transport can 
lead to chemical abundances that depart from local ther- 
modynamic equilibrium. This phe nomenon has long, 
been established in Jovian planets l |Bgjghay_fc Lewis I 
19781 iFeglev fc Prinn I 119851 iFeglev fc Lodders I 119941 
Lodd ers fc Feglev 111994ft . In brown dwarfs, the strongest 
expected signatures of departures from equilibrium 
chemistry are an overabundance of CO and depletion 
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of NH3. The pr ediction of non-equ ilibrium chemistry 
in brown dwarfs f Feglev & Loddersl 11996ft was quickly 
confirmed by the detect ion of the 4.7 /im band of CO 
in Gl 229B ijNoll et al. Ill997ft . and subsequently by the 
detection of NH3 in a high-resolution near-IR spectrum 
ijSaumon et alJl2000D . Nonequilibrium chemistry may 
be common in cool brown dwarfs since M' photometry 
shows a systematic flux depletion in T6 - T9 dwarfs that 
can be attributed to en hanced CO 4.7 /im absorption 
( Goli mowski et, al"1l2004fl relative to models in chemical 

equilibrium that 

With a spectral type of T7.5 ijBurgasser et al. 112006ft . 
Gl 570D is one of the latest T dwarfs known and its IRS 
spectrum shows the strongest 10-11 ixm NH3 absorp- 
tion observed so far ijCushing et al. ll2006ft . It is a par- 
ticularly interesting target since it is a wide companion 
to the well-studied K4 V star Gl 570A and to a pair of M 
dwarfs, Gl 570BC. The distance, metallicity, and age of 
the system very effectively reduce the parameter space 
allowed to model Gl 570D. iGeballe et a!~l IplOl used 
these constraints to determine the effective temperature 
T e ff and gravity g of Gl 570D and to analyze its optical 
and near infrared spectrum. Here, we take advantage of 
the new spectroscopic information provided by the IRS 
spectrum to redetermine (T e g,g). We also perform the 
first quantitative analysis of the NH3 10 - 11 /jm band 
for any brown dwarf. 

2. OBSERVATIONS AND DATA REDUCTION 

Our analysis of Gl 570D is based on spectroscopic data 
covering most of its spectral energy distribution. We 
combine a Keck LRIS 0.63 - 1.01 /an optical spectrum 
l|Burgasser et al. 1 1200^) with a UKIRT CGS4 near IR 
spectrum cover i ng 79 - 1.35 /im and 1.43 - 2.52 ^m 
ijGeballe et al. 1 12001^ . and a IRS Sho rt Low module 
spect rum which covers 5.43 - 14.68 /xm (|Cushing et al. I 
2006). Absolute flux calibration for the UKIRT and 
IRS spectra was obt ained by using MKO photometry 
l)Leggett et al. 112002ft and IRAC band 4 photometry (B. 
Patten, priv. comm.), respectively. The far red optical 
spectrum was then scaled to match the flux density level 
of the flux-calibrated UKIRT spectrum. 
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3. MODELS: EVOLUTION AND SPECTRA 

We analyze the spectrum of Gl 570D with a combina- 
tion of model atmospheres, synthetic spectra and evo- 
lution sequences. Its T7.5 spectral type strongly sug- 
gests that its atmosphere is free of condensates and we 
use cloudless model atmospheres and sy nthetic spectra 
which were most recently described in iFortnev et al. I 
(2006) and references therein. Aver aging recent metal- 
licity determinations for Gl 570A dThore n fc Feltzing I 
l2000t ISalitos et al. 11200,4 IValenti fc Fisher 1 120051) gives 
[Fe/H]=0.09±0.04. This slight enrichment is negligible 
in view of the other uncertainties in the data and models 
and our entire analysis assumes solar metallicity. 

We have developed an evolution code for brown dwarfs 
that assumes an adiabatic i nternal structur e . It u ses the 
H/He equation of state of iSaumon et al. I (|2005). Nu- 
clear ener gy generation through th e first branch of the 
pp-chain (Burr ows fc Liebert 1119931) is included. We use 
the NACRE nuclear reaction rates (jAngulo et al~lll999|) 
and ion-ion and electron-ion screening factors (Chabrier, 
priv. comm.) The surface boundary condition is ob- 
tained from our atmosphere models. All brown dwarf 
atmosphere models become convective at depth, where 
the specific entropy of the gas becomes essentially con- 
stant. The value of the entropy at the bottom of the 
atmospheres S a t m gives the entropy for the adiabatic in- 
terior. The boundary condition is described as a tab- 
ulated function S'atmCToff , <?, [Fe/H]). Thus, our cooling 
seq uences are fully consistent w ith the atmosphere mod- 
els (|Chabrier fc Baraffe 111997}) . Details of the evolution 
calculation will be provided in a future publication. 

4. DETERMINATION OF THE PHYSICAL PARAMETERS 



We follow the method developed bv ISaumon et al 



( 2000) and initially applied to Gl 570D bv lGeballe et al. 



(2001). Briefly, we used the optical, near-IR and mid- 
IR spectra to determine the integrated flux observed at 
Earth (1.875±0.107x 10" 12 erg/s/cm 2 ). 9 This represents 
about 70% of the bolometric flux of Gl 570D. We use our 
grid of model spectra to compute the bolometric correc- 
tion and the bolometric luminosity, L^ ol (T c g , g), which 
depends on the yet- to-be-determined T c g and g. An in- 
dependent relation L£ ol (T e ff , g) is given by our evolution 
sequences. A consistent solution between the spectrum 
and the evolution is obtained by imposing L| ol = L^ oV 
This gives a family of solutions T e s(g) (Fig. 1). A com- 
bination of age indicators s hows that the Gl 570 system 
is between 2 and 5 Gyr old (jGeballe et al. ll2001j) , which 
constrains the range of the gravity to that shown by the 
isochrones in Fig. 1 and the models in Table 1. Mod- 
els A and C bracket the age limits and model B cor- 
responds to the midpoint. This new determination of 
the physical paramet ers of Gl 570D differs from that of 
iGeballe et al. I ()200lD in several ways. The observed SED 
has been greatly extended by the Spitzer IRS spectrum 
and the near-IR spectrum has been recalibrated using 
the latest MKO photometry. The atmosphere models 

9 The uncertainty estimate is obtained from the uncertainties 
in the absolute calibration of each piece of the SED (3% for the 
optical spectrum, 5% for the near-IR spectrum and 6.5% for the 
IRS spectrum). A total uncertainty of ±5.7% is based on conser- 
vatively assuming that the individual calibration uncertainties all 
add in the same direc tion. The 1% uncertainty on the distance 
IPerrvman et al. 1997) is negligible by comparison. 



have steadily improved, most notably with the inclusion 
of a more extensive line list for C H4 i|Champion et al.1 
H992tlFreedman fc Lodders l2006|) . Finally, the evolution 
is computed with an independent code, using our cloud- 
less atmosphere grid as the surface boun dary condition, 
rather than relying on the sequences of iBurrows et all 
()1997l) . Despite all these changes, the parameters ob- 
tained for Gl 5 70D a re remarkably similar to those of 
IGeballe et all l|2001fl (crosses on Fig. 1). Most of the 
difference can be attributed to a shift of the isochrones in 
the (Tpf f , g) plane between ou r evolution sequences and 
those of IBurrows et alJ l)1997[) caused mainly by the dif- 
ferent surface boundary conditions. This shows that the 
method is quite robust and that further improvements in 
the models will not affect the values in Table 1 signifi- 
cantly. 

Figure 1 shows that the main uncertainty in the 
(T e g,g) of Gl 570D arises from the uncertainly in the 
age of the system. By comparison, the uncertainty on 
Lboi from the flux calibration of the data is small. 

Given T c g, g, the radius R(T c g,g) from the evolution, 
and the parallax of Gl 570A, we can compute model 
fluxes at Earth that can be compared directly to the 
data. The model B spectrum is shown in Fig. 2. This is 
not a fit of the spectrum per se, but it is obtained solely 
from the procedure outlined above. The agreement with 
the optical and near-IR data is generally excellent. The 
deviations in the 1.6 - 1.7/xm region are due to the in- 
complete CH4 line list. The reason for the underesti- 
mated K band flux is unclear. Models computed with 
the slightly higher metallicity of Gl 570A would give a 
higher K band flux. A similar effect could also be ob- 
tained by flux redistribution from a stronger CH4 1.6 /im 
band computed from a more complete line list. We find 
that a cloudy model with a thin cloud dec k parameterized 
by a sedimentation parameter f s Br i = 4 l|Fortnev et al. I 
2006; lAckerman fc Marlev l l2001). can also improve sub- 
stantially the agreement with the observed K band peak. 

In the mid-IR, the model systematically underesti- 
mates the flux beyond 9 /im, which coincides with strong 
absorption by NH3. The outlier models (A & C in Ta- 
ble 1) give very similar synthetic spectra. The differ- 
ences remain below ±0.03 mJy in the near-IR and below 
±0.2 mJy in the mid-IR, with higher gravity models giv- 
ing lower fluxes. 

5. NONEQUILIBRIUM CHEMISTRY OF NH 3 

Vertical transport in the atmosphere of a T dwarf will 
result in a reduced NH3 abundance, which will increase 
the model flux for A ij 9 /im and improve the agreement 
with the IRS spectrum. We consider atmosphere models 
that depart from chemical equilibrium to optimize the fit 
to the data. 

The kinetics of the chemistry of carbon and nitrogen 
in giant pla nets and brown dwarfs has been described 
in det ail s in iLewis fc Priim" il980t) : Feelev fc Lodders 
< 199 1 : IGriffith fc Yelle I (119991): ILodders fc Feglev 
1 2002j) : iBezard et al. I (|200^) : lLodders fc Feglev I i|200fl) 
and the effect on th e spectra of brown dwarfs in 
ISaumon et al~l l)2003|) . We only provide a brief 
summary of nonequilibrium chemistry here. 

5.1. Nonequilibrium carbon and nitrogen chemistry 
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The chemistry of carbon and nitrogen in brown dwarfs 
atmospheres is essentially described by the net reactions 

CO + 3H 2 <-> CH4 + H2O (1) 

N 2 + 3H 2 ^2NH 3 . (2) 

Because of the large binding energies of CO and N 2 , both 
reactions proceed much more slowly to the right than to 
the left. Vertical mixing dredges up hot gas that is rela- 
tively rich in CO and N 2 to the cooler part of the atmo- 
sphere, where CO and N 2 are converted very slowly into 
CH4 and NH3, respectively. The net result is a relative 
overabundance of CO and N 2 in the upper atmosphere 
and decreased abundances of CH4 and NH 3 in the upper 
atmosphere. The abundance of H 2 is also reduced by 
conservation of the total number of oxygen atoms. 

The correct chemical pathway and reaction time scale 
for the conversion of CO into CFL remains uncertain 
(lYung et al. Il988t iGriffith fc Yellell999HBezard et al I 
l2002HLodders fe Fedev ll2002HVisscher fe Feelev II2005L 
We ad opt the "fast" chemical time scale of lYung et al. I 
(1988). Fo r the N 2 convers i on, w e adopt the time scale 
given in Lod ders fc Feglevl l|2002j) . 

5.2. Vertical Transport 

Slow mixing can occur in the radiative zones of atmo- 
spheres through a variety of processes such as the tur- 
bulent decay of waves propagated upward from the con- 
vective/radiative boundary and instabilities arising from 
rapid rotation. In the absence of any detailed modeling 
of these processes, we adopt a simple model of verti - 
cal mixing by eddy turbulence ijGriffith fc Yelle I 1X9991. 
This process is analogous to diffusion and occurs over a 
characteristic time scale r mix ~ H 2 /K, where H is the 
pressure scale height and K is the coefficient of eddy dif- 
fusion. K ranges from ~ 10 2 to 10 5 cm 2 /s in planetary 
stratospheres and it is the only free parameter in our 
non-equilibrium models. 

The time scale for mixing in the atmospheric convec- 
tion zone is r conv ~ H c /v c where H c and v c are the 
convective mixing length and velocity, respectively, and 
are evaluated with the mixing length theory. For conve- 
nience, we choose H c = H . Note that r conv << r mix . 

5.3. Quenching of the Chemistry 

Th e chemical abundance profiles ijLodders fc FegleTI 
2002) in model atmosphere B are shown in Fig. 3. Depth 
in the atmosphere is indicated by the local temperature. 
The equilibrium abundances are shown with dashed lines 
that overlap solid lines for H 2 and CH4. The transi- 
tion in the carbon chemistry from CO deep in the atmo- 
sphere to CH4 near the surface occurs at logT ~ 3.25. 
The N 2 to NH3 transition occurs around logT ~ 2.9. 
Remarkably, and in contrast to the carbon chemistry, 
the N 2 /NH3 equilibrium ratio remains nearly constant 
deep in the atmosphere due to a near cancellation of the 
increased temperature (favoring N 2 ) and the increased 
pressure (favoring NH3) along this portion of the atmo- 
sphere profile. 

The chemical time scales for the conversion of CO 
and N 2 are shown by the heavy solid and dotted lines, 
respectively. These vary over more than 20 orders of 
magnitude throughout the atmosphere and become ex- 
tremely long in the upper atmosphere. Mixing occurs 



rapidly in the convection zone but it is about 7 orders of 
magnitude slower in the radiative zone for our choice of 
K = 100cm 2 /s. In regions where r chcm < r mix , chemical 
equilibrium prevails. Where r chcm > T mix , however, abun- 
dances of CO, CH4, N 2 , NH3 and H 2 will depart from 
equilibrium (solid lines). The nonequilibrium mole frac- 
tions are approximately determined by their values where 
T rhrrn = r mix . I n our calculation, we adopt the scheme of 
Smith (1998) to determine this "quenching level." As 
can be seen in Fig. 3, there can be more than one level 
where r chcm = r mix , resulting in alternating zones in and 
out of chemical equilibrium. For simplicity, we consider 
only the uppermost (lowest T) crossing. It turns out 
that in all cases, this is an excellent approximation to 
the full solution for the chemical profile since the deeper 
quenching occurs well below the level of formation of the 
corresponding molecular bands. 

Figure 3 shows that above the CO quenching level 
(logT ~ 3.03) the nonequilibrium abundance of CO is or- 
ders of magnitude higher than at equilibrium. The abun- 
dances of CH4 and H 2 are barely affected for this choice 
of K (slow mixing). The NH3 abundance is reduced by 
about one order of magnitude in the upper atmosphere. 
In this model, N 2 is quenched below the bottom of the at- 
mosphere model, in a region where N 2 /NH3 is still nearly 
constant. The fact that the equilibrium mole fraction of 
NH3 is nearly constant for logT <; 3.1 implies that the 
non-equilibrium NH3 abundance profile is completely in- 
sensitive to the choice of the eddy mixing coefficient K 
over a very wide range (log K ^ —8). This means that all 
plausible choices of K will give the same NH3 abundance 
and therefore the same mid-IR spectrum. 

5.4. Optimal model 

In view of the above remarks, we arbitrarily choose 
K = 100 cm 2 /s to compute nonequilibrium spectra for 
the three models in Table 1. Besides the choice of equi- 
librium or nonequilibrium chemistry, the only fitting pa- 
rameter left is the gravity, constrained to be between 
the values for models A and C (Table 1). We determine 
the goodness of fit of a synthetic spectrum by computing 
the x 2 between the model and the data for A > 9 /im, 
which is where NH3 features are seen in the spectrum 
(Fig. 2). For this purpose, model spectra are smoothed 
to the IRS Short Low point source spectral resolution of 
A/ A A = 126 with a Gaussian filter and then binned to 
the instrument's wavelength sampling. The uncertainty 
(a) on the fitted \ 2 is obtained by fitting 5000 simulated 
data sets generated by adding a random Gaussian noise 
distribution to each pixel with a width given by the indi- 
vidual pixel noise (~ 0.1 mJy). We also fit renormalized 
spectra to account for the ±6.5% calibration uncertainty, 
which turns out to be the dominant source of uncertainty. 
Nonequilibrium model C always gives the best fit. A rea- 
sonable fit (2.8cr worse than model C) is obtained with 
nonequilibrium model B if the observed flux is increased 
by the full 6.5% uncertainty. The best equilibrium model 
fit is 6.7ct worse than the best nonequilibrium model fit. 

The nonequilibrium model C spectrum is shown by the 
dashed (blue) curve on Fig. 2. The agreement with the 
data is remarkable considering that all we have done is 
include non-equilibrium chemistry to the calculation and 
increase the gravity slightly, with no other free param- 
eter. We conclude that the NH3 abundance is reduced 
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by nonequilibrium chemistry and that the parameters of 
Gl 570D fall between models B and C, and much more 
likely closer to the latter. 

6. ALTERNATIVE EXPLANATIONS FOR THE NH 3 
DEPLETION 

We briefly consider five possible explanations that 
could account for the enhanced mid-IR flux in Gl 570D: 
a low nitrogen abundance, an error in the temperature 
profile of the atmosphere, uncertain NH3 opacity, time 
variability, and photodissociation of NH3. 

6.1. Low nitrogen abundance 

The metallicity of Gl 570D is very nearly solar, based 
on abundance determinations of Na, Si, Ca, Ti, V, Cr, Fe, 
Co, and Ni in Gl 570A that range from [X/H]= -0.10 to 
0.24, with a typi cal spread of ~ 0.1 between independent 
determinations ij SjmjQs e^_al. l l2005l iThoren fc Feltzing I 
l2000HValent,i fc Fisher II2005[K To the best of our knowl- 
edge however, the N abundance has not been measured. 
A low [N/H] value would result in a decreased NH3 abun- 
dance that could potentially explain the relatively high 
mid-IR flux. The dependance of the mole fraction of am- 
monia on the nitrogen abu ndance is log Xnh 3 ~ 0.5[N/H] 
l)Lodders fc Feglevll2002ft . In our models, the strongest 
NH 3 absorption occurs at a depth of logT = 2.65 in the 
atmosphere where the non-equilibrium chemistry reduces 
-X"nh 3 by a factor of 10 (Fig. 3). To obtain a similar re- 
duction from a low nitrogen abundance would require 
[N/H]= —2, which would be an extreme departure from 
the abundances of other metals in Gl 570A. 

In principle, N-bearing condensates may also cause a 
depletion of N-bearing gases in the atmosphere. How- 
ever, there are no N-bearing condensates expected along 
the (T, P) profile of Gl 570D that can significantly re- 
move nitrogen, and temperatures are not low enough to 
allow NH4SH or NH3 condensation as in the much cooler 
upper atmospheres of the Jovian planets. 

6.2. Uncertainties in the temperature profile of the 
atmosphere 

Our determination of Lboi and our solution for (T e g, g) 
are in excellen t agree ment with the values obtained by 
IGeballe et al~l l|200lD without the benefit of the IRS 
spectrum. This indicates that while the near-IR flux has 
remained constant, the unknown mid-IR flux of Gl 570D 
in March 2000 could not be significantly different from 
that measured in Feb 2005. The fact that our equilib- 
rium model is in very good agreement with the full SED 
of Gl 570D (Fig. 2) supports the notion that all pieces 
of the SED are consistent with each other. Finally, T c g 
would have to increase to ~ 1050 K to obtain an equilib- 
rium NH3 mole fraction similar to the non-equilibrium 
value needed to fit the IRS spectrum. Such a high T e g is 
completely incompatible with all of the Gl 570D observa- 
tions. If we consider the temperature profile of the atmo- 
sphere, an increase of over 200 K (~ 40%) at P 1 bar 
would be required; an implausible error in our models, 
an inconsitent with the near-IR spectrum of Gl 570D. 

6.3. Uncertainties in the NH3 opacity 

The opacity of NH3 is compute d from a line list ob- 
tained by combining the GEISA ijHusson et al. I Il992]l 



and HITRAN l|Rothman et al. 112004}) data bases, com- 
plemented with recent labo ratory measurements and 
theoretical calculations (see IBurrows et al.l lll997|) for 
more details). To our knowledge, this compilation of 
NH3 opacity is the most complete presently available. 
Nonetheless, the line list is only complete up to T = 
300 K and transitions arising from levels that are ex- 
cited only at higher temperatures are underrepresented 
or missing. Qualitatively, the NH3 opacity we calculate 
for Gl 570D is a lower limit to the actual opacity, imply- 
ing that the derived underabundance of NH3 is an up- 
per limit. More quantitatively, the 9-14 /xm band of NH3 
arises from the ground state (i/q — ► V2), with a small over- 
laping contribution from the v-i excited state (yi — > 2^2) . 
The energy differences between other excited vibrational 
states of NH 3 is too large for the corresponding hot bands 
to overlap with the Vq —> v% band. We are confident that 
the line list for the 9-14 /im band on which our analysis 
of the abundance of NH 3 is based is quite complete. 

6.4. Time variability of Gl 570D 

The spectroscopic data analyzed here was gathered 
over a period of several years: The optical spectrum is 
from observations take n in March 2000 and Feb 2001 
llBurgass er et al. 1l2003l) . t he near-IR spectrum in March 
2000 l|Geballe et al. I200ll) and the Svitz er IRS spectrum 
in February 2005 l|Cushing et al. 1 12006^> . It is conceiv- 
able that the luminosity of Gl 570D may have varied 
during this 5 year period, affecting our determination 
of T e g and g and fortuitously result in our excellent fit 
of the combined spec trum with a model depa rting from 
chemical equilibrium. iBurgasser et al. 1 l)2003j) do not re- 
port any variation in the optical flux of Gl 570D be- 
tween their two observations taken one year apart, and 
the reported equivalent width of Cs I and Rb I lines 
show no variation within 1.5<r. The near IR spectrum 
is cal ibrated to JHK MK O photometry taken in Feb 
2000 ijGeballe et al. Il200lj) . Repeat observations in July 
2005 are fully consistent (< la) with the original pho- 
tometry. The 2MASS photometry of Gl 570D was ob- 
ta ined in May 1998. Applying the color transformations 
of iStenhens fc Leggett I l|2004l) shows agreement to 2a 
(J), la (H) and < la (K) with the MKO magnitudes. 
Based primarily on the stability of the near-IR MKO 
photometry, we conclude that any time variability in Gl 
570D over a period of 7 years is at most at the 2-3% level. 
This is less than our estimate of a 5.7% uncertainty on 

Nevertheless, we consider the possibility that the abso- 
lute mid-IR flux may not be representative of the average 
or typical state of Gl 570D. If we relax the absolute flux 
calibration of the IRS spectrum, do we still see a NH3 
depletion? For this purpose, we have repeated the fitting 
procedure described in §5.4, with the flux calibration ad- 
justed by a factor that minimizes the value of x 2 - Allow- 
ing for the estimated noise in each pixel of the observed 
spectrum, we find that the non-equilibrium models al- 
ways fit better than the equilibrium models at the ~ 17<r 
level, with insignificant variation between models A, B 
and C. 

6.5. Photodissociation of NH3 

NH3 is a molecule that is easily dissociated by UV 
photons ~ 1900 A. Since Gl 570D is a companion to a 
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main sequence star, it is possible that the UV flux from 
the primary results in a depletion of N Hs in Gl 570D. 
This p rocess was discussed in details by iSaumon et al.l 
( 2000) in the case of another late T dwarf with a com- 
panion, Gl 229B, where it was found to be completely 
negligible at the level of the photosphere. By comparison, 
the UV flux of the dK4 star Gl 570A is about 10 times 
larger than that of the dMl star Gl 229A (bas ed on the 
NextGen models of lAllard k. Hauschildtl lj!999|) 'l but the 
flux incident on Gl 570D is smaller bec ause of the much 
great er projected separatio n (1525 AU llBurgasser et al. I 
12000ft rather than 44 AU l|Golimowski et al. 111998(1 for 
Gl 229B). Since the atmospheres of both brown dwarfs 
are quite similar, we can expect that the photodissocia- 
tion of NH 3 in Gl 570D will be 10 x (44/1525) 2 ~ 0.008 
times that in Gl 229B. The Gl 570 BC pair of M dwarfs 
lies about 10 times closer to Gl 570A than Gl 570D and 
their contribution to the UV flux incident on Gl 570D is 
a small correction to the above estimate. 

7. CONCLUSION 

Our method to determine T c g and the gravity of brown 
dwarfs is most effective when applied to objects with 
known parallax and which are members of systems where 
the primary provides the metallicity and age indicators. 
The method is based on integrated fluxes and evolution 
sequences rather than a direct fit of the spectrum, an 
approach that largely eliminates biases due to remaining 
systematic problems in the model spectra. With recal- 
ibrated optical and near-IR spectra and, most impor- 
tantly, a Spitzer IRS spectrum, we have redetermined 
(T R ff, g) for Gl 570D to find nearly the same values as 
lOehalle et aFI ll200l . 

An excellent fit of the IRS spectrum can only be ob- 
tained if the NH3 abundance reduced by nearly one or- 
der of magnitude in the upper atmosphere. This is 
most readily explained by vertical mixing that drives 



the nitrogen chemistry aw ay from chemical equilibrium 
l)Lodders fc Fegfe v"l 120021k We have considered several 
other possible explanations but none of them are plausi- 
ble. 

We find that in late T dwarfs, the noncquilibrium 
N2 /NH3 ratio generally does not depend on the efficiency 
of mixing in the atmosphere. Therefore, while the 10 - 
11 /im band of NH3 is a sensitive indicator of noncqui- 
librium chemistry in brown dwarf atmospheres, it pro- 
vides little information about the efficiency of mixing. 
On the other hand, the CO abundance is very sensitive 
to the choice of the mixing parameter K, and photomet- 
ric or spectroscopic observations of the 4.7 ^tm CO band 
is the most effective way to determine K . Our study pro- 
vides the most accurate parameters of any known brown 
dwarf (Fig. 1): T cS = 800 - 820 K, log g = 5.09 - 5.23, 
R/Rq = 0.0901 - 0.0853, log L / Lq = -5.525 to -5.551, 
M/Mj = 38 - 47 and an age of 3 to 5Gyr. The upper 
limit of this range (Fig. 1) is at least 2.8cr more likely 
than the lower limit. Gl 570D joins Gl 229B as the second 
T dwarf where the effect of nonequilibrium chemistry due 
to vertical transport has been demonstrated explicitly. 

This work is based (in part) on observations made with 
the Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under NASA contract 1407. T. Roellig and M. Mar- 
ley acknowledge the support of the NASA Office of Space 
Sciences. K. Lodders is supported by NSF grant AST04- 
06963. Part of this work was supported by the United 
States Department of Energy under contract W-7405- 
ENG-36, and NASA through the Spitzer Space Telescope 
Fellowship Program, through a contract issued by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with NASA. This research has ben- 
efited from the M, L, and T dwarf compendium housed 
at DwarfArchives.org and maintained by C. R. Gclino, 
J. D. Kirkpatrick, and A. J. Burgasser. 
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Fig. 1. — T e f{ and gravity for Gl 570D. Our brown dwarf evolution tracks are shown by heavy black lines labeled with the mass in Mq. 
Isochrones (blue dotted lines) are labeled in Gyr. The nearly vertical lines (solid and dashed, red) are the locus of (T e g ,g) points with 
losL/L(n = — 5 .525 d z 0.024. The solid dots show the three model of Table 1. For comparison, the crosses show the values obtained by 
Gcballe ct al. (2001). [See the electronic edition of the Journal for a color version of this figure.] 
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Fig, 2 . — Comparison of the obse rved and modeled SEP of Gl 570D. The data (heavy black curve and histogram) are from Bureasscr et al. 
120031) . lG"eballe et alH 120011) and ICushing et al. I 120061) . The equilibrium model B (Table 1) is shown by the thin (red) solid line, and 
the best fitting nonequilibrium model (model C with log K = 2) is shown by the dashed (blue) line. The main molecular absorbers are 
indicated. The models are plotted at a spectral resolution of R=500 (upper panel) and R,=100 (lower panel), approximating the resolution 
of the data. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 3. — Chemical profile of model atmosphere B (Table 1). Mole fractions of H2O, CH4, N2, NH3 and CO are shown in equilibrium 
(dashed curves) and out of equilibrium (log .ft" = 2, solid curves) as a function of the temperature in the model atmosphere. The pressure 
increases from the top of the atmosphere (left) toward the bottom (right). Heavy black lines show the mixing time scale {r mim ) and the 
time scale for the destruction of CO (tco, solid) and N2 (tn 2 , dotted). Time scales are given on the right axis. The mixing time scale is 
nearly discontinuous where the atmosphere becomes convective (logT ^ 3.2). The 10—11 (im NH3 band is formed in the logT = 2.65 to 
2.9 region. [See the electronic edition of the Journal for a color version of this figure] 
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TABLE 1 

Range of physical parameters of Gliese 570D 



Model 






log L/Lq 


Mass 


Radius 


Age 




(K) 


(cm/s 2 ) 


(Mj) 


(Re) 


(Gyr) 


A 


782 


4.95 


-5.503 


31 


0.0950 


2.0 


B 


800 


5.09 


-5.525 


38 


0.0901 


3.2 


C 


821 


5.23 


-5.551 


47 


0.0853 


5.0 



